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We suggest that a previously conjectured relation between Structure Functions (SF) for 
nuclei and nucleons also links distribution functions (df) for partons in a nucleus and in 
nucleons. The above suggestion ensures in principle identical results for SF F^, whether 
computed with hadronic or partonic degrees of freedom. In practice there are differences, due 
to different F% input. We show that the thus defined nuclear parton distribution functions 
(pdf) respect standard sumrules. In addition we numerically compare some moments of nu- 
clear SF, and find agreement between results, using hadronic and partonic descriptions. We 
present computations of EMC ratios for both and compare those with hadronic predictions 
and data. In spite of substantial differences in the participating SF, the two representations 
produce approximately the same EMC ratios. The apparent correlation between the above 
deviations is ascribed to a sumrule for F%- We conclude with a discussion of alternative 
approaches to nuclear pdf. 
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^ ! Abstract 
P 

We suggest that a previously conjectured relation between Structure Functions (SF) for nuclei 
and nucleons also links distribution functions (df) for partons in a nucleus and in nucleons. The 
above suggestion ensures in principle identical results for SF F^, whether computed with hadronic 
or partonic degrees of freedom. In practice there are differences, due to different F^ input. We 
show that the thus defined nuclear parton distribution functions (pdf) respect standard sumrules. 
In addition we numerically compare some moments of nuclear SF, and find agreement between 
results, using hadronic and partonic descriptions. We present computations of EMC ratios for 
both and compare those with hadronic predictions and data. In spite of substantial differences in 
the participating SF, the two representations produce approximately the same EMC ratios. The 
apparent correlation between the above deviations is ascribed to a sumrule for F^- We conclude 
with a discussion of alternative approaches to nuclear pdf. 
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I. INTRODUCTION. 



A large body of data is presently available on inclusive lepton scattering from nuclei and 
on subsequently extracted nuclear Structure Functions (SF). Standard treatments of those 
data use hadron degrees of freedom, whereas in an alternative approach, one employs df of 
partons in nuclei. In one of the first papers on EMC ratios, Akulinichev et al related df for 
quarks in a nucleon, nucleons in nuclei and quarks in nuclei by a generalized convolution, 
schematically written as [l| (see also Refs. 

fq/A = fN/A* fq/N, (1-1) 

Variations of the above and various approximations have since been the preferred tool for 
the analysis of EMC ratios, whether in terms of the hadron or a parton representation. 

There exist alternative treatments of nuclear pdf. For instance Frankfurt and Strikman 
discuss nuclear pdf, but do not relate those to the pdf of a nucleon jj . Other approaches 



parametrize information on SF at some sea 



e Qq, and then evolve to a desired Q 2 . Data 



then determine the introduced parameters p,|6|]. Next we mention quark models for nuclei, 
which have been used in direct calculations of df. Their potential is presumably limited to 
the lightest nuclei []]]. Finally there are approaches, where the effect of a nuclear medium 
on a nucleon or a quark is replaced by mean fields p, 0, [llj . 

In the present note we make a simple, nearly natural choice for nuclear pdf, which is 
free of adjustable parameters. Those satisfy basic sum rules and produce the same as 
computed from a hadronic base. In the conclusion we compare some of the above-mentioned 
alternative proposals with our choice. 

II. A FEW ESSENTAILS. 

We start with the cross section per nucleon for the scattering of unpolarized electrons 
with energy E over an angle 9 

d2 ° A ^J V) = °m(E; 6, v) [^F^x, Q 2 ) + Q 2 )tan 2 (£/2)] (2.1) 

om is the Mott cross section and F-^ 2 (x, Q 2 ) are the standard nuclear SFs per nucleon. Those 
depend on the squared 4-momentum transfer q 2 = —Q 2 = — (\q\ 2 — u 2 ) and on the Bjorken 
variable < x = Q 2 /2Mv < A in terms of the nucleon mass M. 



Next we make explicit the specific relation (jl.lj) between nuclear and nucleonic SF 
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F%' n are SFs of p, n, whereas F% defines a nucleon SF, obtained by weighting F%' n with 
Z,N. 5N/A denotes the relative neutron excess. 

Our approach draws on the Gersch- Rodriguez-Smith (GRS) theory for inclusive scattering 
of non-relativistic projectiles [3], in which case the linking f PN > A is the SF of a fictitious 
nucleus, composed of point-like nucleons. Eq. (j2.2j) results from a covariant generalization 
of the GRS theory 13]. The theory is non-perturbative with on-mass shell nucleon SFs. 

Eq. (jl.lj) has originally been formulated in the Bjorken limit We have postulated (j2.2j) 
;o hold for finite Q 2 |ll| and an extensive body of data in the r ang es x > 0.2 ; Q 2 > 2.5 GeV 2 



14, 



15J] appears to be accounted for by the relation (f2~2|) 

m 



Eq. (j2.2j) does not contain a term, which account for quarks in virtual mesons 19( , and it 
also lacks (anti-)screening effects [^J. Those are negligible for x > 0.2 and we limit ourselves 



to that range 



44|. 



In the following it appears useful to separate the nucleon SF F, 



N _ rpN,NE , W N,NI 



into nucleon elastic and inelastic components, which correspond to absorption processes 
of a virtual photon on a N, 7* + iV — > N, (NE), or 7* + iV — > (hadrons,partons) (NI). 
Elastic components for a N are proportional to the standard combinations of squared electro- 
magnetic form factors G% M (Q 2 ) and vanish, unless x = 1. With [G^] 2 = [Z{G P ) 2 + 
N(G n f}/A, one has 
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The corresponding nuclear NE (QE) components are from Eq. (j2.2j) seen to be 

f PN > A (x,Q 2 ) 
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(2.4) 
(2.5) 



(2.6) 
(2.7) 
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In particular for the lightest nuclei the normalized f PN < A peaks around x « 1, and the 
same holds for the above QE components F k ' NE (x, Q 2 ). The above summarizes elements of 
a hadronic description of nuclear SF: we now turn to a partonic representation. 



III. A SIMPLE CHOICE FOR NUCLEAR PARTON DISTRIBUTION FUNC- 
TIONS. 

We start with the leading order twist contributions to the dominant NI components of 
nucleon SF for finite Q 2 . With no danger of confusion we shall omit one or both arguments 
x,Q 2 . Neglecting heavy quark contributions and decomposing quark df q = q v + q into 
valence and sea quarks parts, one has 



Fl = | (Au v + d v + 8u + 2d + 2s' 
F 2 n = | (u v + Ad v + 2u + 8d + 2s) (3.1) 



Similarly for a 'nucleon', defined as the Z, N weighted p, n, one has 



= dixqi = 777 u v + d v + 2u + 2d + - ^^-{u v - d v + 2u - 2d) , (3.2) 
*— ' 18 L 5 oA J 

Nuclear pdf ought to reproduce nuclear SF Fj^, just as proton pdf do for F%, Eq. (|3.1|) . 
but that necessary requirement is insufficient for their determination. In a construction we 
proceed in two steps. First we choose F^(q A ) to be the same combination of df of partons 
in a nucleus, as the above F 2 , Eq. ()3.2|) . is for the averaged nucleon 



4 35N 



F 2 A = Yaaqt = u A + d A + 2u A + 2d A + V - ~zr~r~ (u A - d A + 2u A - 2d A ) (3.3) 
L — 4 18 L 5 i>A J 

i 

Upon substitution into Eq. ([2.2)1 . and using Eq. ()3.2|) . one finds 



4 A 35N 
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Comparison with Eq. ()3.3|) still does not result in a unique expressions for each individual 
pdf q A . Then guided by Eq. (|2.2j) . we make the following second choice, with no flavor 
mixing of valence and sea quarks 

xqf(x,Q 2 ) = j\zf PN >\z,Q 2 )(^u v (- z ,Q 2 



xqf (x,Q 2 ) = jy Z f PN > A (z,Q*)(^)d v (- z ,Q 
xg A (x,Q 2 ) ee [ A dzf PN ' A (z,Q*)(*)g(*,Q 2 

J x \Z / \ Z 



(3.5) 



Eqs. (|3.5|) with one f for all partons (in case for the PWIA) had already be suggested by 
Berger and Coester 2|. 

In view of the meager experimental information on non-valence parton distributions in 



nuclei, we shal 
medium (cf 



also investigate changes when non-valence df are not affected by the nuclear 
), thus 

q A ee q; s A = s A = s; g A = g N , (3.6) 



Eqs. ()3.3|) . ()3.5|) (or alternatively Eq. ()3.6|) ) manifestly produce the same F A in the parton 
and the hadronic representation, provided one uses exactly the same input f PN > A and F% 
in both. 

In practice this is not the case, in particular not for F^ in F^\ In the absence of 
direct information, the CteQ parametrizations exploit data on F®. Using the 'primitive' 
approximation F^ix) = 2F®(x) — F%, F^x) replaces _F 2 D as input. However, the above 
approximation deteriorates with x, increasing from x > (0.25 — 0.30). As we shall shortly 
demonstrate, its use for larger x leads to misfits with data. 

A 8 to the option, in the version CteQo Q we selected the one, with * ciosest to 
the Arneodo parametrization of resonance-averaged data 23|. This appears possible for 
x ^ 0.6 — 0.7. Using SU3 symmetry one obtains in the usual way F£ in terms of g^. 

In contrast, in the hadronic approach one stays as close as possible to data. For not 
too high Q 2 , one can exploit parametrizations of actual F% data (2^]. However, the range 
Q 2 > 3.5 GeV 2 of our interest borders the limits of validity of the data parametrizations, 
we chose to use the above-mentioned Arneodo representation of i 7 ^ j^]- As to F^-, it has 
been obtained for allx by an indirect extraction method (see for instance Ref. 2^)- Fig. 1 
for C = F^/F^ in the two representations, Q 2 = 3.5, 5.0 GeV 2 , makes manifest the above- 
mentioned differences, which increase with x (4^. 

A second subtle difference between the representations is the validity of Eq. ()2.2j) and 
(13. 5|) . The former one is explicitly limited to the nucleonic part of F A , i.e. roughly for 
x > 0.2, while the latter is conceivably valid out to lower x (cf. Ref. 
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The above described freedom in the choice of df for partons in nuclei is not unlimited, 
since df are constrained by sum rules. We first check those, which are directly related to the 
normalization of f PN ' A ; 

a) For any linear combination C of df for valence quarks 



a r i 

A _ A^nN 



dxC A = / dxC N (3.7) 
o Jo 



/ dx(u A + d A ) = / dx(u v + d v ) = 3 
Jo Jo 



(3.8) 



= f ix ^±jk = l (3. 9 ) 
'o 3 J 3 

Eqs. (|3~7|) . (jSSJ) are for C A = u A + d A = u A -u A + d A -d A , respectively C A = (2u A -d A )/3. 

Those express the conservation (per nucleon) of the number of valence quarks in nuclei 

(baryon number) and of charge. 

b) For any linear combination C = a u u v + otdd v of nuclear pdf 

C A (0,Q 2 ) = C N (0,Q 2 ) 
C A (x ,Q 2 ) « C N (x ,Q 2 ) (3.10) 

By construction, Eq. (|2.2jl holds only for the contribution of partons from nucleons, and 
somehow the same is the case for (|3.5|) . For those parts one proves the above equality for 
x = 0, whereas for 0.18 > xq > 0, Eq. ()3.10|) is an approximation. 

c) For any combination xC(x) = x ^ aiqi(x, ) or x J2i a i<li( x ) 

[ A dxxC A (x)= [ A dzzf PN > A (z) f dttC N (t) (3.11) 
Jo Jo Jo 

[ dxx[u A + d A + 2u A + 2d A + 2s A + g A } x = [ dzzf PN ' A (z) * 
Jo Jo 

/ dtt [u v + d v + 2u + 2d + 2s + g] t 
Jo 



I dtt [u v + d v + 2u + 2d + 2s + g] t (3.12) 
Jo 



Again, Eq. (J3.12)) is a special case, related to the momentum sumrule, which does not 
exactly carry over to the nuclear case. However, the peaking of the normalized f PN > A causes 
the ^-integral in Eq. (|3.11j) to be very close to 1. In the case of the nuclear momentum 



sumrule Eq. (J3.12j) . the deviations are really minute (see Ref. |28| for an entirely different 
way to mend the momentum sumrule violation in the PWIA). 



The incompleteness of the above-mentioned Eq. ()2.2j) for x < 0.2 does not in practice 
constitute a problem. For one, the conservation of the number of valence quarks per nucleon 
is guaranteed by unitarity, i.e. the normalization of f PN ' A . As the momentum sumrule 
(J3.12j) illustrates, more is requires than unitarity, in case the fact that the widths Ax(Q 2 ) 
of the peaked / are appreciably less than the support of x. Comparing the two sides of the 
expressions in Eq. ()3.12|) for 4 targets and Q 2 = 3.5, 5.0 GeV 2 , we find differences of no more 
than « 1%. 

The above sum rules involve pdf in a nucleus and in a nucleon. Various sumrules involve 
moments of nuclear SF and one may compare numerical results in the hadronic representa- 



tion (see for instance Ref. 



15 



29 



30]) and in the partonic one. For instance, Eq. (|3.3|) and 
(I3.12|) hint, that there exists a sumrule for M\ = F A /x = Y2 a iQ^- The contributions of 



valence quarks is 5/6, but all other parts diverge (see for instance 31, 32). Those divergences 
cancel in differences of any pair of those ratios. Thus from Eqs. ()3.3|) . (|3.6|) (for simplicity 
we disregard 5N/A corrections) 
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(3.13) 



dx(u A + d A ) - 

'0 JO 

For all A' ytz A (including A' = 1 with a x-support, maximally different from (0,A) for 
a target A), the upper integration limits in Eq. (|3.13|) are unequal. However, for the 
above-mentioned reason, one may neglect the contributions to the integrals in ()3.13|) for 
x < 0.20 and x > 0.95. Hence there effectively is an approximate common upper limit 
x\j ~ 1 <C A, A', beyond which the difference (F A — F A ') is negligible and the same holds 
for xl < 0.18. For the pair D, Fe and the chosen three Q 2 values, the above difference of 
the integrals in the hadron and pdf representations is ~ —0.03 [46]. 

Of special interest is the zeroth moment of nuclear SF, which is related to the momentum 
sumrule. Using Eq. ()2.2|) 



Mi 



dxFt 



A-NI , 
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(3.14) 



The same moment of iso-singlet NI parts in the pdf representation is from Eqs. ()3.3j) . ()3.6|) 
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and (|3.8|) seen to be 



M A = / dxF A ' NI (x) = 4 / d.r.r 



1 

5 

5 
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u A + d A + 2(u A + d A ) + ts A 



J dxx{^ 



6 



I - / dxx[-s A + g J 



(3.15) 



1- / cte(-s(x) = / rfxF 2 (7V) (x) (3.16) 



NE parts are small for the considered Q 2 . When included, the normalization of f PN < A 
guarantees M A,NE = Mq' ne , i.e. NE parts of M A are also A-independent. 

In Figs. 2a-2e we show differences of valence, sea quark and gluon distributions functions 
in a nucleus and for the p. We choose 5 targets and display results only for Q 2 = 5 GeV 2 , 
because for the Q 2 range considered there is hardly any Q 2 -dependence. Differences increase 
with increasing mass number, and change sign at roughly x = 0.2 and 0.8. 

To the extent that nuclear sea and gluon distributions are close to the nucleonic ones, Eq. 
(I3.16|) shows that in the pdf representation, M A is practically A- independent. Emphasis is 
on the standard ~ 50% reduction of the nucleon valence contributions 5/18 due to gluons 
(see for instance |3JJ), which carries over to nuclear df. 

It is not feasible to directly verify the hadronic result ()3.14j) . which requires F A to be 
known over the entire relevant x-range. Present data are for x > 0.3 for the lowest Q 2 (for 
which Eq. (J3.14J) is not accurate) and x > (0.5 — 0.6) for medium Q 2 . For both the missing 
information for x < 0.5 contains the major contribution to M A . We thus take recourse to 
nuclear SF, computed from Eq. (j2.2j) . Wherever data are available for the DIS x-range, 
those appear to agree very well with computed ones. 

Table I displays the zeroth moments M A , A = {D, He, C, Fe}, computed in the hadronic 
and the partonic representations (Eq. ()3.14j) . respectively ()3.15j) . ()3.16|) ). Comparison of 
columns 2 and 3 clearly shows that the shape f PN ' A effectively cuts the upper limit of 
long-range integrals at x\j ~ 1.0. 

The lowest moments M A {Q 2 ) are seen to be dependent on A in a weak and not smooth 
fashion. Going from D to Fe, those moments are 3-4, respectively 8% smaller than the same 
for the averaged nucleon Mq N \q 2 ): D and He clearly do not follow the smooth behavior of 
all other nuclei. One also notices that Mq are weakly descending functions of Q 2 . 

At this point we return to a verification of the result ()3.14|) . In spite of the lacking data 
on F A for small-x, it is possible to reach an indirect, approximate experimental verification 



()3.14|) in the hadronic representation. The method exploits knowledge on F A for small x and 
the smoothness of the same in the region x > x m , x m being the above-mentioned, smallest 
x measured [3, I35I. One then interpolates nuclear SFs in the intermediate region, where 
data are missing |25J] and subsequently approximate calculates the lowest moments M A . In 
spite of the fact, that the 'missing' region contributes the major part of the integrals (|3.14j) . 
the moments turn out to be surprisi ngly close to the ones computed above. Forthcoming 
data from JLab experiment E03-103 |36( will enable a sharpening of the above method. 

The approximate A-independence of the above zeroth moments of any pair of SF entails 
the same for the differences 

M A - M A ' = [ A dxF 2 A -' NI - [ A dxFf> NI « f U dx[F^ NI - Ff' m ] « (3.17) 

JO JO Jx 

The vanishing of the difference is attributed to about equal effective x-ranges (xq,xu) ~ 
(0.15,0.85), which replace actual unequal supports j^j]. Eq. ()3.17j) implies that in the 
above common interval, the difference of two SF has to change sign at least once, or in 
different terms: the generalized EMC ratios fj, A ' A '(x, Q 2 ) = F A {x, Q 2 ) / F A ' (x, Q 2 ) passes 
the value 1 in the above x-interval, as is indeed observed for all A, Q 2 (see Ref.|38j for a 
discussion in an entirely different context). 

The simplest cause for approximate A- independence of Mq would be the same for F A , 
but that appears not to be the case. In the dominant classical region x < 0.90, differences 
F A — F A ' in both representations grow with x beyond « 0.18 and may become as large as 
50-60 %, which is far larger than the spread in M A (cf. Table I). We shall shortly return to 
this point. 



IV. EMC RATIOS IN THE PARTON DISTRIBUTION FUNCTIONS REPRE- 
SENTATION. 



We have computed F A ' NI in the pdf representations, using Eq. ()3.3j) . ()3.5j) and (JUS}. To 
those we added the NE components (|2.7j) . which are only relevant for |1 — x\ < 0.05. The 
total EMC ratios F A I FP are then compared with recently determined counterparts in the 
hadron representation [38]. For the range 0.2 < x < 1.2 those, and avaibale data for He, Fe 



and Au 



39 



40 



41[, are shown in Figs. 4a,b for Q 2 = 3.5, 5.0 GeV 2 



Up to x < 0.65 there is close agreement between the computed ratios. Beyond that 
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point deviations set in, which grow while pdf ^> hadr. The latter situation is reversed for 
x > 0.90. Both representations overestimate the relative maxima in fi A around x = 0.9, but 
over the entire x-range, the hadronic results are closer to the data than those for the pdf. 
We attribute this to inferior F^ input for larger x, which propagates into F% through the 
use of Eq. (|2.2j) . Also of interest are the slightly lower pdf results for x < 0.65 and the much 
higher ones beyond up to x ~ 0.9. 

The apparent insensitivity of EMC ratios fi A to the representation, in spite of the large 
differences in the participating F A , combined with points a), b) in Section III suggest the 
following: Irrespective of the cause of the dependence of nuclear SF on A and/or repre- 
sentation, the approximate independence of the zeroth moments (13.14)) on both, forces the 
differences in F A in the regions x < 0.18 and 0.18 < x < 0.90 to be nearly balanced. Con- 
sequently, if EMC ratios in one area are in some way ordered in A, that ordering is inverted 
in the second one. It almost seems that deviations of EMC ratios from 1 can be generated 
by an integral-preserving, afhne transformation with x ~ 0.18,0,90 as fixed-points, and 
having in particular for A < 12, a characteristic A-dependence, and having x ~ 0.18,0,90 
as fixed-points. 

Finally, out of sheer curiosity we followed the pdf results down to x = 10~ 5 . Results for 
F A hardly change from their ^-independent values around x = 0.15, causing EMC ratios to 
stay close to 1 for decreasing x. The above is actually observed down to x ~ 10~ 3 . Only for 
the smallest x do screening effects deplete df and cause /i A to slowly reach values ~ 0.6 — 0.7. 

V. COMPARISON, CONCLUSIONS. 

In this note we defined df for partons in a nucleus, using for the latter a relation between 
nucleonic and nuclear hadron SF without any adjustable parameter. We showed that either 
exactly or very closely, those respect basic sumrules. Observables such as F A may be 
expressed in either representation, and by construction produce in principle the same F A 
for identical input F% ,n . 

We have shown that this is not the case in practice, specifically not for F£. The pdf choice 
rests on the 'primitive' approximation F% = 2F 2 D — F$, which is satisfactory for x < 0.3, 
but deteriorates with increasing x, while a well-founded extraction method has been used 
in the hadronic representation. Consequently EMC ratios, computed in the two versions, 
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practically coincide for x < 0.65. For large x deviations appear, which reflect the similarly 
deviating F% from a better founded extracted function: the hadronic representation of EMC 
ratios produces the better fits to the data. 

The above, and the fact that deviations of EMC ratios from 1 appear in distinct areas to 
be balanced for all A, seems linked to the lowest moment M A of F A : those are for any Q 2 
practically independent of A. 

We continued pdf calculations down to the smallest x, which is the region were the above 
criticism does not hold. However, there Eq. ()3.5|) misses primarily (anti-)screening effects. 
Nevertheless, down to x ~ 10~ 3 the agreement with data persists, but the pdf results cannot 
describe anti-screening depletion of df in \i A for the smallest measured x. 

Our almost natural choice of distribution functions of partons in a nucleus is clearly one 
out of many possible ones, and we shall mention a few suggested alternatives. For instance 
Eskola etal address participating nuclear SF in EMC ratios, which are generated from 
parametrized input for a reference Q%, and which are subsequently evolved to the desired 



- H. ' 



Q |5j . Parameters are constrained, for instance by fixing the average position of minima of 
EMC ratios. Data determine those parameters. 

Next we mention Kumano and co-workers who, in spite of different x-support for F A 
and F N , assume a linear relation between df for partons in a nucleus and in a nucleon 



.1 



The species-dependent, relating weight functions Wi(x, Z, A) contain parameters for 
a scale Qq and the resulting df of partons in a nucleus are again evolved to any desired 
Q 2 , ultimately producing parametrized nuclear SF and EMC ratios. A large number of 
adjustable parameters leads to fits from the smallest x up to x < 1. With the connecting 
weight function having no meaning beyond x — 1, the interesting region x > 1 is out of 
reach in that approach. A more serious drawback of the method may be the lack of physical 
meaning of the weight functions Wi and its parameters. 

Finally we discuss approaches, where df of partons in a nucleus are those for a nuc eon 
bound in scalar and vector mean fields [si, 0], which couple to quarks in a nucleon 42^ . 
Off-hand, the above and our phenomenological approach seem to have little in common, but 
this is actually not the case, and it is instructive to trace the connection. We recall that 
the original proposal Eq. ()2.2|) was inspired by a model, where valence quarks in a nucleus 



cluster in bags 



llj . Total interactions between quarks in two different bags have there been 



replaced by phenomenological NN forces, acting on the centers of those bags, thus replacing 
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quark dynamics by those for hadrons. The above treatment of a nucleus partially calls back 
quark degrees of freedom, and are from Eq. (|3.4j) seen to mix the latter with hadronic 
degrees of freedom through the SF f PN ' A . 

In contrast, the hybrid meson-quark coupling model, where the interactions of a single 
quark in a given nucleon with all quarks in the remaining A — 1 bags are replaced by mean 
fields. As expected from its intermediate position, it is possible to sum in that model the 
above meson-quark interactions over the valence quarks in a bag, and to construct a NN 
interaction with many-body components, mediated by the same mean boson fields. This 
has recently been shown to be possible 1431 . 

nn 

The above has been applied to nuclear matter [8|, |9( and it is clearly of interest to see 
applications of finite nuclei. A first example is a treatment of 3 He in the PWIA [10] and one 
should look forward to results for higher A. 
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